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Abstract 

An analysis is presented that examines the modulation of different 
instability modes satisfying the triad resonance condition in time and space 
in a three dimensional boundary-layer flow. Detuning parameters are used for 
the wavenumbers and the frequencies. The nonparallelism of the mean flow is 
taken into account in the analysis. At the leading-edge region of an infinite 
swept wing, different resonant triads are investigated that are comprised of 
traveling crossflow, stationary crossflow, vertical vorticity, and Tollmien- 
Schlichting modes. The spatial evolution of the resonating triad components 
are studied. 
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I. Introduction 


An important stage in the transition from laminar to turbulent flow is a 
region of nonlinear development, characterized by a broad spectrum of non- 
linear ly interacting disturbances. The character of the nonlinear development 
is strongly dependent on the initial spectrum of the disturbances. 

Most theoretical and experimental work was performed using the Blasius 
boundary layer to investigate the instability mechanism that break down into 
turbulence. On the other hand, little is known about the physical phenomena 
that leads to transition in cases like swept wings where the boundary layer is 
three-dimensional (3D). In this situation, the boundary-layer profile 
consists of a streamwise velocity component in the direction of the external 
inviscid flow and a crossflow velocity component normal to it along the wing 
surface. Due to that, different types of instability modes may exist and 
different possible interactions may occur between these modes resulting in 
stability characteristics that is much different from what linear theory would 
suggest. 

The resonant interaction of three waves is considered one of the mecha- 
nisms that play an important role in determining the nonlinear characteristics 
of the development of disturbances, leading to transition. Such resonance 
occurs whenever the real wavenumbers k and frequencies 0 ) satisfy the 
conditions 

k i ± k£ ± k3 = 0, o)| ± 0)2 i 0)3 = 0 
with corresponding signs being taken. 
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In Blasius boundary layer, there usually exist triads comprised of 

two-dimensional wave propagating in the flow direction and two-obliquely 

12 3 

propagating plane waves. Raetz ' and Stuart established the occurence of 

triad resonances for certain neutrally linear stable waves. Craik ~ examined 

the occurence of this triad resonance for certain unstable waves over a flat 
8 

plate. Lekoudis derived the nonlinear spatial and temporal evolution 
equations of the triad waves by relaxing Craiks assumption of perfect 
resonance. The resonance model of wave interaction is one of the 

9 

experimentally observed phenomena (e.g. r Kachanov et al. , Kachanov and 
Levechenko 10 , Saric and Thomas 1 X , Saric et al. 12 ), in two-dimensional flows 
that play an important role in the 3D secondary instability that break down 
into turbulence . 

However, in 3D boundary layers, as on a swept wing, the triad is 

comprised of three resonantly interacting 3D waves that may propagate in 

different directions. Because 3D boundary layers are usually rich in 

instability modes, one expects the possible evolution of different triads with 

i 3 

different interacting modes that resonate. Lekoudis confirmed the existence 
of a triad on a swept wing that consists of three unsteady crossflow modes, 
but the interaction coefficients and the amplitudes of the interacting waves 
were never calculated in 3D boundary-layer flows. 

It is known that transition prediction methods used for modern LFC trans- 
port depend primarily on the use of the e 11 criterion. This method is based 
only on the exponential growth of small individual disturbances within a 
boundary layer according to linear stability theory. With the interaction of 
the linear modes of these disturbances and the possibility of rapid growth of 
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their amplitudes, as is the case for Blasius boundary layers 14 ' 15 , the e n 
criterion is no longer valid and a modified one is needed. 

In this article, we investigate the evolution of resonant triads in 3D 
boundary layers. The triads investigated are comprised of different 
instability modes, stationary crossflow (CF), traveling crossflow, vertical 
vorticity (W) , and Tollmien-Schlishting (TS) modes. In section II the 
nonlinear analysis of the triad resonant interaction is developed. The mean 
flow used in the calculations is the boundary layer on a modern LFC transonic 
23° swept wing. In our analysis the growth of the boundary layer is taken 

39 

into account assuming that it is of the same order as the nonlinear effects. 
Details of the mean flow are given in section III. Section IV discusses the 
numerical procedures. Results for different mode-mode interactions are given 
in section V for parallel flows, while nonparallel flow effects are discussed 
in section VI. Then we end with concluding remarks. 

II. Nonlinear Analysis 

We consider the nonlinear interaction of wave packets in a 3D incompres- 
sible boundary layer on a swept wing. The flow field is governed by the non- 
dimensional incompressible Navier-Stokes equations. The Cartesian coordinate 
system used has the x-axis in the direction of the normal chord, the z-axis 

along the span, and the y-axis normal to the surface. The Reynolds number 

i/2 

R = u*L*/v* is based on a reference length L* = (v*s*/U*) t where s* is 
e e e e 

is the distance along the airfoil surface surface, and v* is the kinematic 
viscosity coefficient evaluated at the edge of the boundary layer. 

The mean flow is assumed to be slightly nonparallel, with e a small 
parameter characterizing the flow divergence, and identified with 1/R. The 
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method of multiple scales 16 is used to introduce the slow scales (X|,Zj.,t|) = 
(ex,ez,et) that govern the growth of the boundary layer, the modulation of the 
disturbance amplitude, and the change in the eigenfunction. While the phase 
of the disturbance changes over the scales x, z, and t. 

To determine the wave packet solution of the governing equations, we 
assume that the flow quantities possess uniformly valid expansions in the 
form. 


* r n 3 

u =* U(x 1 #y,z <| ) + A e ^ n (x,x 1 ,y,z,z 1 , t,^ ) + 0(e ) 

I n=1 

2 n 3 

v = ev(x .y,z ) + I e v (x,x ,y,z,z ,t,t ) + 0( e ) 

n I 11 

n=1 

n 3 

w = W(x ,y,z ) + l e w (x,x ,y,z,z ,t,t ) + 0( e ) 

II n 1 11 


p = P(x ,z ) + l e n p n (x,x ,y,z,z ,t,t ) + 0( e ) 
n=1 


(1) 


( 2 ) 


(3) 


(4) 


Where U, V, W, and P are the steady mean-flow quantities and u, v, w, and 
p are the unsteady small disturbances . To account for the simultaneous 
effects of the flow divergence and the disturbance nonlinearity, we assume 
that both effects can be expressed by the same expansion parameter e. 

Substituting equations (1) - (4) into the governing Navier-Stokes equa- 
tions after transforming the time and space derivatives, subtracting the mean- 
flow terms, and equating the coefficients of like powers of e, we obtain. 
Order £: 


3u 3v 9w^ 

VWV = aT- + ~ = 0 


3u 3u 


3z 

3U 


3U 1 3P 1 -1 2 
VWV = TT + 3^1 + W 3i" + 3^“ R V U 1 = 0 


(5) 

( 6 ) 


4 


(7) 


3 1 1 


3v 1 . , 


£l + 

3p i 

i p“ 


»- 

4J 

\<ro 

>zr + 1 

9z 

3y 

T XV 

v Vi 

3w i 

3w i 

9w 


, 3W 1 

3p 1 

’l } = 3t“ 

+ U 3x + Sy^ 1 

+ i 

*3z~ 

+ 3z 


- R -1 7 2 w. 


( 8 ) 


where 


v 2 = ^i- + ^- + i- 

3x 2 3y 2 3z 2 


Order e 


9u , 5w , 

WW = - 3^-- 3^ 


(9) 


2 

3u 3u 3u 3p 3 u 

L 2 (u 2 ,v 2 ,p 2 ) = - 3 ^ - - ” 3 ^ - 3 ^- + 2R ' < 


9 2 u, 


3x3x l 3z3z^ 


3u 3u 3u 3U 9U 1 3U , 

- (u i 3 T +v i aT +w i aT 1 ’ 3x i 1 3y 3z i 1 


(10) 


3v 1 

L 3 (V 2' P 2 ) = " 3T 


8v 


9v 


3 2 v 


3 2 v 

j 1 W 1 . OR " 1 ( L_ + 

^ ^7 3x3 Xl 3z3z 1 


1 


) 


, 3W 1 3W 1 3v i, 3v i 3V . 

- (U 1 JT + V 1 IT + w i TT ] - [v ay- + 37 V 1 ] 


( 11 ) 


VW'V 


3w i 3w i 


3t 


1 3X 1 
9v* 


9u i _ ^ a2 "i 

Vz~ + 9x9x^ + 3z3z^ 


1 3w i 3w i 3w 

* (U 1 3T + V 1 3T + W 1 3^’ " 1 3^1 " 'W - ^ a*,"! 


3w i 3W 

+ V-~- + ^-*,1 (12) 


Here, the leading order problem governs the linear wave in a parallel 
flow, while the higher order problem includes both nonparallel and nonlinear 
effects . 


5 



2a. First-order equations 


Consider the nonlinear interaction that may exist among three wave 
packets centered at the frequencies o^, 0 ) 2 , and 0 ) 3 . Thus we express the 
solution of equations (5) - ( 8 ) as a linear combination of three interacting 
waves according to, 

3 ie 

q 1. * J WW t »n < V y ' z t ,e " <”> 


where q 1 m ,m= 1,...6 stands for u± , Suj/ 3y, v L ,w L , 3w x /3y, and p A respec- 
tively and 

30 

rr m WV (,4) 

ae 

7z~ m e n (x r z i ) (15) 

ae 

j~ = -u n , n=1 ,2,3 (16) 


The phase functions 6 n are assumed to be continuously differentiable, that 
is 


3a 


3z. 


36 


3x_ 


(17) 


The a and 0 are the complex wave numbers in the x and z directions 
n n 

given by a = a + ia . and 0 =0 + i 0 . , and to are the complex 

n nr ni n nr ni n 

frequencies given by a) = o) + io) . . They satisfy the resonance conditions, 

n nr ni 


o) — 0) - a) = ea 

3 r 1 r 2 r t 


a, - a -a = ea 
3 r 1 r 2 r x 


0 - 0 - 0 = ec 

3r P 1r P 2r z 


(18) 

(19) 

( 20 ) 
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where the detuning parameters a , a and a [all 0(1)] are introduced to 

t x z 

introduced to express quantitatively the nearness of the above resonance. 

Substituting (13) -(16) into equations (5)-(8), separating coefficients 

i6 n 

for e , n=1,2,3, and writing the result as six first-order systems of 
ordinary differential equations, we obtain. 


DC - l (b ) C = 0, m=1 , . . .6 
mn j — m 3 n 3 n 


subject to the boundary conditions 


'in ' c 3n - * 0 “ y * 0 


'in' Sn, C 4n * ° 


as y + » 


(21) 


(22) 

(23) 


where D = d/dy, and the nonzero coefficients of (b . ) are given in 

mj n 

Appendix A. 

2b. Second-order equations 

In order to determine the conditions for the elimination of secular terms 
in the second-order equations (9)— (12) , and hence determine the amplitudes 
a n in (13), we seek a particular solution for the second-order equations in 


| the form, 

i 


2m 


3 i9 n 

i ,W 1 "VV t i )e + c - c - 

n=1 


(24) 


where q , m=1,...6, stands for u ,3u /3y,v ,w , 3w /3y and p , respec- 
2m 2 2 2 2 2 2 

tively. Substituting (13)-(20) and (24) into (9)— (12) and separating the 

i 0 

coefficients of e n , n-1,2,3, we obtain three separate systems of equations, 
where each can be written as six first-order sets of equations in the form, 


DiJ> 


mn 


6 

- I 

j=i 


(b .) 

m] n jn 


I , m= 1 , ... 6 
mn 


(25) 


i 
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subject to the boundary conditions 


* 


In 

In 


= + 


3n 


•*3n'V 


Mn 
+ 0 


0 at y = 0 
as y + 00 


( 26 ) 

(27) 


The coefficients (b .) are the same as in equation (21), while the 

m} n 

inhomogeneous terms I are functions of the first-order eigensolutions £ 

mn mn 

of equation (13), a , 6 , w and the mean-flow quantities. The terms I 

n n n mn 

include all nonparallel and nonlinear effects. 

Since the homogeneous parts of equation (25) are the same as in equation 
(21), and since the latter has a nontrivial solution, the inhomogeneous equa- 
tion (25) has a solution if, and only if, and only if, the inhomogeneous parts 
are orthogonal to every solution of the corresponding adjoint homogeneous 
problem; that is, 

» 6 

/ l * S* dy = 0, n=1,2,3 (28) 

' „ mn- mn 

0 m=1 


where are solutions of the adjoint problems, they are, 

n 


6 


D?* - l (b . ) 
mn . ** mi n 

3 = 1 

C* 

3n 

— 0 , rri— 1 , ... 6 

(29) 

? 2n = ? 4n = ? 5n = 

0 

rt 

u 

o 

(30) 

5 2n' 5 4n' C 5n * ° 


as y + 00 

(31 ) 


for n= 1 ,2,3, and b . = -b. . 

mj j m 

The solvability conditions (28) give the following differential equations 
for the evolution of the amplitudes a^ in time and space. 


3a 3a ■ 3a^ __ 

h 11 at 1 + h 21 13T + h 31 37" + £ 1 h 4i a i + £ 2 h 5i a 2 a 3 e = ° 


(32) 
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9a 9a 9a — r 

h 12 3t + h 22 3x” + h 32 3z + e i h 42 a 2 + e 2 h 52 a i a 3 S 


( 33 ) 


9a 3 3a 3 9a 3 _r 3 

h . _ TT-- — + h . _ -r — + _ *r + £ h ^ _ a _ + £_h_ _a^ a _ e 

13 3t 23 3x 33 9z 1 43 3 2 53 1 2 


= 0 


(34) 


where h„ , h n , h_ , h. , h c , and T are given in Appendix B, and ( ) 
in 2n 3n 4n 5n n 

indicates a complex conjugate of ( ). 

Equations (32) -(34) account for the combined effect of the growth of the 

boundary layer, and the nonlinear interactions. The parameters e L and e 2 

are shown here ( £| = e 2 = e in the analysis) only to indicate terms due to 

each effect. If e 2 << e if the nonlinear interaction can be neglected. When 

e l << e 2' the nonparallel effect can be neglected. It is worth noting that 

with this formulation both nonlinear and nonparallel effects are independent 

of the particular disturbance quantity being considered. 

Solution of equations (32) -(34) for a general initial condition is not a 

simple task. However, by assuming the spatial modulation of a single 

frequency disturbance on an infinite span wing, a situation closer to 
1 7 i 8 

experiments r , we can allow for modulation only in the x-direction, and 
equations (32) -(34) can be simplified by 
9 3 


9t 9z 


= 0 


U) — (U — ur 
3r 1 r 2r 


S 3r " P 1r ” P 2r 


0, a t = 0 


o, a =0 
z 


(35) 

(36) 

(37) 


It is convenient to introduce the transformation 


A = £a exp (-fa.dx-B.z) 
n n ‘ nx m 


(38) 


then equations ( 32) — ( 34) reduce to 


dA _ 

ST" (G 1 - a ii )A 1 + H 1 A 2V 


(39) 
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( 40 ) 


id) 

ST’ la 2 - a 2i )h 2 * VlV 


= (G - a . )A + H, A As" 1 * 
c 3 3i 3 3 12 


where 


G n = -e^, (42) 

H - - r-— -, (43) 

n n 

2n 

4> * J ea x dx ( 44) 

In equations (38), 0^ is a parameter and is made equal to zero in equations 

( 39 ) — ( 41 ) • We note that for the case of an infinite span wing, if the initial 

wave has 0^ = 0, it will remain zero downstream. 

In order that the interaction coefficients and the wave amplitudes 

are uniquely defined, it is necessary to specify the normalization imposed 

on the eigensolutions £ of the first-order problem. This is chosen such 

mn 

that the maximum of the r.m.s. of r (c = u ) over y is equal to one. 

In In n 

Note that the nonparallel coefficients do not depend on the normalization 

of C or • 

mn mn 

To derive real equations for the amplitudes and phases, we let, 

i * 

A - -k e n (45) 

n 2 n 


* L X r 


G = G e 
n n 


H = H e , n=1 ,2,3 
n n 


Substituting (45) -(47) into (39)- (41) and dropping the *, we obtain. 


dT = (G 1 COS X 1 “ °1i )A 1 + lV2 A 3 COS (Y + T 1 + ♦ ) 
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(49) 


1 

dT = (G 2 COS X 2 ' a 2i )A 2 + 2Vl A 3 C ° S (Y + X 2 + ♦ ) 


dA 


*T = (G 3 COS X 3 “ Q 3i )A 3 + lVl A 2 C ° S ( ' Y + T 3 ' ♦ ) 


(50) 


dy 


dj = (G 3 sin x 3 - 0 2 sin x 2 - G, sin X, ) * t « 3 ^ . . . . 


A 1 A 2 

sin (-y + T, - <(> ) 


A A A A 

- H 2 -LI sin (Y + T 2 + ♦) - H 1 sin (y + T, + + )] 


(51) 


where 


Y " S * X 2 - X 1 


(52) 


III* Mean Flow 

The mean flow used in these calculations is the boundary layer with 

suction on a 23° swept infinite span wing. The airfoil section (designated 

SCLFC ( 1 ) -051 3F ) is supercritical with normal chord c = 6.44 ft. This wing 

was the subject of extensive experiments designed to examine supercritical 

laminar flow control technology at the Langley Research Center ' • Linear 

2 1 22 

stability calculations for this wing have been given by El-Hady ' , 

2 3 24 

Mack , and Berry et al. 

Frees tream conditions for the present calculations are Mach number = 0.82 
and a chord Reynolds number of 20 x 10 2 * * * 6 . The upper surface pressure coeffi- 
cient distribution is shown in Fig. (1) together with the suction parameter 
distribution, and the distribution of the boundary layer maximum crossflow 
component |v n | maximum along the chord. 

The three-dimensional boundary-layer solution is calculated using a 

25 

boundary layer program that is adapted from the code of Kaups and Cebeci for 
laminar, compressible boundary layers with adiabatic and wall suction boundary 
conditions. The code assumes zero pressure gradient along the wing generator. 
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IV* Numerical Procedures 


For n » 1, 2, and 3, the set of equations (21 )— (23) and their adjoints 
( 29)-( 31 ) can be solved analytically in the freestream at y = y , producing 
three linearly independent, exponentially decaying solutions with the 
characteristic values. 


A. = -(c£ + eV /2 


(53) 


-[a 2 + 8 2 + iR(a + W 8 - u )] 1/2 

n n n e n n 


A. = A. 


(54) 

(55) 


with the freestream solution as initial condition, equation (21) is integrated 

26 

from y - y to y = 0 at the wall, using a variable step-size algorithm , 
e 

based on the Runge-Kutta-Fehlburg fifth-order forumlas. The solution is 

orthonormalized at a preselected set of points using a modified Gram-Schmidt 

procedure* A Newton-Raphson technique is used to iterate on the eigenvalue to 

satisfy the last wall boundary condition. The eigensolutions associated with 

the adjoint problem can be determined by integrating equations (29)-(31) using 

the same procedures and the same previously determined eigenvalues. 

To evaluate the nonparallel terms 3? /3x„ and da /dx„, for n ** 1, 2, 

mn 1 n i 

and 3, we differentiate the first-order problem (21) with respect to x x and 
obtain. 


D( 


3 ? 6 

mn.) _ l 


3x 


1 


j-1 


3C,_ 


6 

- I 

3-i 


3(b .) 
m] n 

3x. 


V 


m 


35 1n 

3 Sn 

3C 

9x i 

9X 1 

3x 

9? 1n 

35 3n 

9C 4n 

9X 1 # 

9X 1 ' 

9X 1 


at y = 0 


as y + » 


» • • 6 


(56) 


(57) 


(58) 
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The homogeneous parts of equation (56) have a nontrivial solution* Their 
eigenvalues and adjoint are the same as those for the first-order problem 
(21). Then, by applying the solvability condition, we obtain. 


da 

n 




1, 2, 3 


( 59 ) 


where h are given in quadratures in terms of £ , £* , a , and 6 (see 

6n mn mn n n 

Appendix C) . Condition (59) permits the integration of equations (56)-(58) to 

determine using the same procedures as those for the first-order 

problem, but for nonhomogeneous sets of equations* 

With all terms in h and h known, the nonparallel coefficients G 

2n 4n n 

are calculated from equation (42), separately for each wave (n = 1, 2, and 3). 

Whereas the interaction coefficients H are calculated from equation (43) 

n 

using the parallel results for the three waves altogether. 

The calculations are repeated at different streamwise locations to 

evaluate G and H for given waves of fixed physical frequency f in Hz, 
n n 

and fixed physical spanwise wavelength X (nondimensionalized with respect to 

z 

the normal chord c) that satisfy the conditions (36) and (37). For different 
initial amplitudes of the respective waves, the amplitude evolution equations 

(48) -(51) are then integrated using a fourth-order variable interval Runge- 

27 

Kutta method by Fehlberg . 


V. Nonlinear Effects in Parallel Flows 

When the amplitude of the disturbance, although small, is sufficiently 
large, nonparallel effects are thought to have no substantial influence. In 
this section, we assume that e i << e 2 such that nonparallel effects are 
neglected. 
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5a • Interaction between traveling CF modes 


First , we study a triad that is comprised of traveling CF waves near the 

leading edge of the swept wing. The triad has the frequencies f x = 100 Hz, 

f 2 = 200 Hz, and f 3 = 300 Hz. Hie corresponding spanwise wavelengths are 

X = 0.0006, 0.0006, and 0.0003, respectively. Calculations start at R = 260 
z 

(0.078% chord) where the three waves are unstable in the absence of the 
interaction. The frequencies and spanwise wavelengths satisfy the resonant 
conditions (36) and (37) at all streamwise locations. 

When the initial amplitudes of the triad waves are comparable in 
magnitude, the interaction between the waves is found to be very weak (not 
shown). It is also found that the initial spectrums of the triad amplitudes 
and phases play an important role in the interaction process. For example, 
Fig. (2) shows the modulation with R of both A 3 and the phase angle y 
defined by equation (52). Hie initial amplitude A 30 = 0.0001, and the 
initial phase angle Yo = 0 (all calculations are for y 0 = 0 unless other- 
wise stated). For (A l0 = A 20 ) < 0.0001, the amplitude A 3 is hardly affected 
by the interaction. With the increase of A^q an d A 2 q, a sharp increase in 
A 3 occurs starting at an earlier streamwise location. Fig. (3) shows the 
modulation of A 3 for the same conditions as in Fig. (2) but for Yo = In 
both cases the modulation of A A and A 2 is almost unaffected by the inter- 
action. Fig. (4) shows that the modulation of A 3 is affected by its initial 
amplitude. If A 30 is large, A 3 is hardly affected by the interaction, 
until later downstream after the amplitudes A^ and A 2 become large enough. 
While for small A 30 , the effect of the interaction shows up early upstream 
with a very sharp increase. 

A strong interaction may also occur and amplify A^ or A 2 depending on 
the appropriate initial amplitude and phase spectrum. Fig. (5) shows strong 
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modulation of A x when Aj^q is very small compared to A 20 and A 3 Q. The 
same picture is almost true for the modulation of A 2 when a 2 q is very 
small compared to A^g and In both situations A 3 undergoes weak 

variations from its linear modulation (not shown). Fig, ( 6 ) gives the depend- 
ence of A | on the initial amplitude A a q when A 20 - 0,0005, and A 3 q = 
0 . 002 . 

Fig. (7) shows the variation with R of the interaction coefficients 
H n' n = 2, and 3, given by equation (43). It shows that | H n |' correspond- 

ing to f 3 = 300 Hz and X = 0.0003, is much higher than those for and f 2 * 

Fig. ( 8 ) shows the variation of the detuning parameter ea x with R. It 
indicates that perfect tuning occurs only at one location. In spite of that, 
preceding results show that a strong interaction continues to exist even if 
the resonant conditions are not tuned provided that the initial amplitudes 
have the appropriate spectrum. 

As we see* from previous calculations, the triad used for this study 

exhibits a strong resonance that may amplify a superharmonic (A 3 ), or may 

amplify a subharmonic (Aj^ or A 2 ). This is again illustrated in Fig. (9). At 

a sufficient distance downstream, the domination of one or the other will 

depend on the spectrum of the initial amplitudes of the interacting waves. 

These results may explain the anomalies found in the crossflow observations of 
28 

Saric and Yeates • In spite of the 1 -cm space streaks they visually 
observed, their hot wire measurements in the boundary layer showed a super- 
harmonic of 0.5-cm periodicity dominated disturbance growth. Reed explained 
this anomalies using an approach that considered the growth of the super- 
harmonic as a secondary instability in the presence of finite amplitude 
unsteady crossflow disturbance. The triad resonance model presented here 
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predicts that a superharmonic (A = 0,0003) will amplify and its amplitude 

z 

ratio will eventually reach three times or more the amplitude of other compo- 
nents of the triad (A « 0.0006/ A = 0.0006) when the initial amplitude of 

z z 

the superharmonic is very small. The same model also predicts that a subhar- 
monic (A = 0.0006) may dominate disturbance growth if its intial amplitude is 
z 

small compared to other components of the triad. 

5b. Interaction between traveling CF and W modes 

In a three-dimensional boundary layer, the disturbance is necessarily 
three dimensional. At the leading edge region of a swept wing and in the 
direction of crossflow, there are really two spectra to be considered to the 
solution of the leading-order equations (21). Crossflow modes (stationary or 
traveling) are given by the eigensolution of equation (21). The same equa- 
tions also admit eigensolutions with v = 0 and p = 0, correspond physically 

to horizontal motions which are called vertical vorticity eigenmodes. These 

3 0 3 3 

eigenmodes, as in the case of two-dimensional flows * , are always damped 

but the damping rate may be quite small so that the nonlinear effect could 

3 2 33 

trigger large instabilities ' 

In our calculations, we were able to converge to a W mode in the 

crossflow direction, that is damped. 

Here, we study the resonant interaction of a triad that is comprised of 

two traveling CF modes (f| = 100 Hz, A^ = 0.0006 and f£ = 200 Hz, A^ = 0.0006) 

and a W mode (fq = 300 Hz, A = 0.0003). Calculations start at R = 779 
0 z 

(2.1% chord), where, in the absence of the interaction, the CF modes are 
unstable while the W mode is damped. Again, the frequencies and spanwise 
wavelengths satisfy the resonant conditions (36) and (37) at all streamwise 
locations . 
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Fig, (10) shows that a W mode with A 30 = 0,0001 resonate with two 
traveling CF modes and becomes strongly unstable in a short distance 
downstream. Fig. (11) shows the influence of A 30 on this instability when 
the CF mode initial amplitudes A^q =* A 20 = 0.002. This strong instability 
may be due to a strong interaction coefficent |h 3 J (corresponding to the W 
mode) compared to | H x | and j 1 r s ee Fig. (12). This strong instability of 
the W mode occurs in spite of the imperfect resonant conditions shown by the 
distribution of the detuning parameter ea in Fig. (13). 

5c. Interaction between traveling CF and stationary CF modes 

Here, we study the resonant interaction that is comprised of two travel- 
ing CF modes (f? = 300 Hz, X = 0.0006 and f 3 = 300 Hz, X = 0.0003) and a 

z ^ z 

stationary CF mode (f x = 0, X = 0.0006) near the leading edge of the swept 

z 

wing. Calculations start at R = 260 (0.078% chord) where the three waves are 
unstable in the absence of the interaction. The frequencies and spanwise 
wavelengths satisfy the resonant conditions (36) and (37) at all streamwise 
locations • 

Fig. (14) shows the effect of the stationary CF vortex with initial 
amplitude A^q = 0.005 on traveling CF modes A 2 and A 3 . Since all waves 
are unstable in the range of Reynolds number considered, the interaction 
process seems to accelerate the growth of A 2 and A 3 compared to their 
linear growth. The detuning parameter distribution is given in Fig. 

(15). Fig. (16) gives the effect of the initial amplitude of the stationary 
vortex on the modulation of A 2 and A 3 • 

5d • Interaction between TS and stationary CF modes 

This type of interaction is a major unanswered question concerning 
swept-wing flows. Some type of interaction between a crossflow vortex and 
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amplifying TS mode is thought to cause premature transition on a swept wing, 

34 

Saric and Reed suggested that the anomaly behavior of transition in the 

3 5 

early LFC work of Bacon et al. when sound is introduced in the presence of 

crossflow vortices can be due to this interaction phenomenon. Also, the 

3 6 37 

experimental work of Poll on yawed cylinders and of Michel et al. on an 

infinite swept wing, observed unsteadiness at transition which might also be 

4 0 

due to this interaction phenomenon. Reed , using a parametric resonance 
model, showed that CF vortices could excite TS modes in the three-dimensional 
boundary layer on the X-21 wing, producing a double exponential growth of the 
TS mode. 

Here, we study one of these possible interactions. In the absence of the 

interaction, a crossflow vortex (f, = 0 and X = 0.003) starts to amplify 

2 

nearly at R = 706 (1.6% chord), then experiences a very slow growth for a 
long distance downstream, until it dies approximately at R = 1663 (13% 
chord). Two TS modes having the same frequency (f 2 = 20 KHz, X^ = 0.012 and 
f 3 = 20 KHz, X^ « 0.0024) will amplify shortly after R = 706 and both decay 
around R = 1150 (5.8% chord). The frequencies and spanwise wavelengths of 
the triad satisfy the resonant conditions (36) and (37) at all streamwise 
locations . 

For different values of the initial amplitude of the crossflow vortex. 
Fig. (17) shows large amplification of the TS modes A 2 and A 3 due to the 
interaction, when A 2 q and A 3 q are small (A 20 - A 3 o = 0.0001). The figure 
also shows a reduction in the vortex amplitude due to the interaction. Higher 
values of the initial TS amplitudes weaken the interaction as shown in Fig. 
(18) for a CF vortex initial amplitude A l q = 0.1. 
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Fig, (19) shows the change in the interaction coefficients with R, where 
j H 3 | and |h 2 | (corresponding to TS modes) are higher than | H i|* Fig- (20) 
shows the distribution of the detuning parameter ea ‘ with R. 

The effect of the initial value of the phase angle y 0 is given in Fig. 
(21) for A10 = 0.1, A 20 = a 3G - 0.0001. The anomaly behavior of A3 at the 
beginning of the interaction is due to different values of Yq . Fig. (22) 
shows how y reaches the same value at some distance downstream in spite of 
different initial values Yo • 

VI. Nonparallel Flow Effects 

When e i = e 2 , nonparallelism of the mean flow comes into play signifi- 
cantly in the nonlinear amplitude modulation. However, the nonparallel terms 
in equations (48)— (51 ) turn out to be most important as the disturbance first 
grows and this, in turn, controls what happens subsequently as the amplitude 
of the disturbance increases. Only the interaction between TS and stationary 
OF modes is investigated in this section. 

Fig. (23) shows the linear parallel and nonparallel amplitude modulations 
of the same instability modes given in section 5d. Fig. (24) shows the 
nonlinear modulation of the amplitudes with the nonparallelism of the mean 
flow included for various initial amplitudes of the interacting modes. This 
figure indicates that during the initial growth or decay of the amplitude, its 
modulation follows exactly the nonparallel development until the amplitudes 
are large enough to interact nonlinear ly. 
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VII, Concluding Remarks 


1 . The preceding calculations show that the development of many triads , 
whose components can, in principle, take part in several resonant interactions 
at once, occurs in three-dimensional flows of boundary-layer type. 

2. An important role in the nonlinear process is played by the initial 
spectrum of amplitude and phases of the triad components. 

3. Due to the interaction of different instability modes, even if they 

are not strongly amplified, the classification concept suggested by 
3 8 

Pfenninger for the stability problem into independent modes is no longer 
valid. 

4. Nonparallel flow effects control the initial development of the 
disturbance triad components while the disturbance amplitude is sufficiently 
small. As the amplitudes increase, nonlinear effects control their subsequent 
spatial development. 

5. The above analysis becomes incorrect with the increase of the 
amplitudes A^, but the nature of a set of phenomena in the 3D boundary 
layers is connected to the resonant mechanism, and may be explained on the 
basis of the results. 
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for n = 1, 2, and 3 
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Fig. 4 The modulation of A 3 with 
R for A x q = 0.0005, A 2 o * 0.002, 
■y 0 = 0, and various values of A 30 
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Fig. 5 The modulation of Aj with Fig. 6 The modulation of A x wit 

R for A x0 = 0.00001, yo a °» and R for A 20 = 0.0005, A 30 = 0.002, 

various values for A 20 and A 30 . y 0 = 0, and various values of A 10 
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Fig. 10 The modulation of a YV mode Fig. 11 The modulation of a VV mode 

with A 30 = 0.0001, yo ■ 0 due to with various Initial amplitude A 30 

resonant Interaction with two traveling due to resonant interaction with two 

CF modes of different Initial amplitudes traveling CF modes A x , A 2 
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Fig. 13 Variation of the detuning 
parameter eo x with R. 
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Fig. 16 Effect of the initial amplitude A i0 of the stationary 
CF vortex on the modulation of traveling CF modes A, and A 3 . 







Fig. 17 The modulation of TS modes A 2 and A 3 due to resonant interaction with 
stationary CF vortex with various initial amplitudes 





Fig. 18 The modulation of TS modes A 2 a 
amplitudes due to resonant interaction witl 
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Fig. 21 Effect of y 0 on the amplification of the TS modes A 2 
and A 3 due to resonant interaction with a stationary crossflow 
vortex A. . 



Fig. 22 The modulation of y with 
R for the conditions of Fig. 22. 
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Fig. 23 The linear parallel and linear nonparallel amplitude modulations of a stationary CF 

mode A x and two TS modes A 2 and A 3 given In section 5d. 
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Fig. 24 The nonlinear nonparallel amplitude modulation of a stationary CF mode A £ and two TS modes 
and Ai given In section 5d for various initial amplitudes of the interacting modes, a) A 10 =.05, A 2 c 
= 0.0001, b) A 10 = 0.01, A 2 q = A 30 = 0.005. Linear parallel ; nonlinear nonparallel . 





